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We calculate the effective interaction W e ff between two holes added to the ground state of 
the repulsive three-band Hubbard model. To make contact with Cooper theory and with earlier 
Hubbard model cluster studies, we first use a perturbative canonical transformation, to generate a 
two-body Hamiltonian. Then, we extend the results to all orders. The approach is exact in principle, 
and we obtain a close analytic expression including explicitly the effects of all virtual transitions 
to 4-body intermediate states. Our scheme naturally lends itself to embody off-site, inter-planar, 
phonon-mediated and other interactions which are not considered in the Hubbard model but may 
well be important. The result depends qualitatively on the symmetry of the two-hole state: 1 £?2 and 
1 Ai pairs are special, because the bare holes do not interact by the on-site repulsion (W=0 pairs). 
The effective interaction in these channels is attractive and leads to a Cooper-like instability of the 
Fermi liquid; however W e ff is repulsive for Triplet pairs. Bound two-hole states of the same nature 
were reported earlier in small cluster calculations by exact diagonalisation methods; only symmetric 
clusters are good models of the plane. Once W e ff is known, the pair eigenfunction is determined by 
an integral equation. We present numerical estimates of the binding energy |A| of the pairs, which 
is in the physically interesting range of tens of meV if unscreened on-site repulsion parameters are 
used. 
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I. INTRODUCTION 

Cooper pairs in high-Tp superconductors jjj are more 
tighly bound and spatially localized than in ordinary 
BCS ones, and this happens in the presence of strong 
correlations. If phonons are mainly responsible for the 
mechanism, an enhanced electron-phonon vertex must 
be involved and a polaron or bipolaron behavior || 
is also a serious possibility. However the wave function 
of the hole system must be such that the repulsion bar- 
rier is minimized and can be overcome by phonon ex- 
change. Several Authors have pushed the argument fur- 
ther, proposing that correlation effects alone could be so 
strong to turn the repulsion into attraction by some sort 
of overscreening and lead to pairing. The proposed elec- 
tronic mechanisms have often been based on Hubbard 
J|,D or t-J models jfj), just to cite a few of the less exotic 
proposals. The phonon-based and electronic mechanisms 
have usually been regarded as alternative, but there is a 
clear need for a theoretical framework where they can 
coexist and can be compared on equal footing; in our 
opinion, all the serious mechanisms must convey part of 
the physics. 

Our first aim here is to propose such a scheme, starting 
with the three-band Hubbard model. The Hamiltonian 
is 

H = H + W (1) 

and the independent hole term reads, in the site repre- 
sentation 

]T [clc d + h.c] (2) 

Cu O n.n. 

where n.n. stands for nearest neighbors. The on-site 
repulsion term will be denoted by 

W = Y^Uin i+ ni-] (3) 

i 

where Ui — U d for a Cu site, Ui — U p for an O. 

Several years ago, pioneers in cluster calculations 
explored the possibility that pairing could result from 
repulsion. They proposed the following definition |7| of 
the energy of the pair: 

A = E(N + 2)+E(N)-2E(N + l), (4) 

where E(N) is the ground state energy of the cluster 
with TV holes, as obtained by exact diagonalisation. A 
negative A means that the ionisation potential decreases 
with increasing the number of holes, and this may be 
taken as an indication of pairing. Several clusters were 
studied for N = 2. It turned out that A < is possible 
only if the off-site U between Cu and O is taken to be 
unphysically large (about 5 eV). Moreover, the physical 
interpretation was not obviuos, and there was no clear- 
cut reason for excluding hole bags. Also, interestingly, 



S. Mazumdar et al. |T^] argued that A < could be an 
artifact due to the neglect of the degrees of freedom of 
the nuclei: when they are left free to move, the state with 
N + 1 holes could gain enough energy from Jahn- Teller 
distortion that A could turn out to be positive, after all. 
Later, we pointed out that a couple of key ingredients 
were missing, namely, symmetry and W=0 pairs. W=0 
pairs are two-hole eigenstates of the kinetic energy Hq 
that are also eigenstates of the on-site repulsion term W 
with eigenvalue 0. There is no on-site repulsion barrier 
to overcome in W=0 pairs. 

The symmetry of the cluster is essential, because only 
clusters with the same point symmetry group C^ v as the 
plane and centered on a Cu ion (fully symmetric clus- 
ters) allow such solutions. Still, also the planar lattice 
structure is essential, because no W=0 pairs occur in 3D 
or in a continuous model, where the interaction term (|^) 
becomes a contact term. We started by the above 
definition of A, and studied the fully symmetric clusters 
with up to 21 atoms by exact diagonalisation; A was 
found to be negative when (and only when) the least 
bound holes formed a W=0 pair. The Hamiltonian had 
been parametrised by electron spectroscopy studies [l^ ] 
and ab-initio calculations [ fL3| , and we used literature 
values. We also considered first-neighbor off-site inter- 
actions The 0-0 hopping can change the order of 
single-hole levels, thereby changing the occupation num- 
ber which is necessary to get a partially occupied degen- 
erate state; however, when the conditions are satisfied, 
the A < behavior results without important modifica- 
tions. The first neighbor off-site Cu-0 interactions were 
found to enhance the effect; the first neighbor off-site O- 
O interactions unpair it, but their effect is small if the 
recommended values of the parameters are used. In the 
present paper, such terms are dropped for the sake of 
simplicity, since they are not essential for our present 
scope. Recent calculations by Schiittler et al |L4|], based 
on a combination of diagrammatic and Quantum Monte 
Carlo methods on an Extended Hubbard Model, support 
the present view that the screened interaction is attrac- 
tive and the attraction is robust against the long-range 
part of the Coulomb repulsion. 

In the fully symmetric Cu-0 clusters a genuine pair- 
ing takes place, due to an effective interaction which is 
attractive for singlets and repulsive for triplets. Let us 
summarize why this conclusion is free from the ambiguity 
pointed out in Ref. fTpf) . 

First, we computed A for our clusters by second-order 
ground state energy diagrams (modified for degenerate 
ground states, when appropriate) [^5| . It is convenient 
to denote the hole orbitals by their symmetry labels, with 
b for b\ and a, (a') for the occupied (empty) orbitals of 
a\ symmetry: we obtained 

A ( 2 ) = _ 2 V W ( a ' h ' x ' x ) 2 _ \^ W(a,a',x,x) 2 

4^ ( £ b ~ £a) (£a' ~ S a ) 

- o a' _ 

(5) 
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where e a is the one-hole energy of the a orbital and so 
on; the sums run over all empty states of the appropriate 
symmetries and involve the matrix elements of the on- 
site interaction W. No contributions arise to second order 
from the empty states of e symmetry since the relevant 
matrix elements vanish. 

Second, we computed to second-order the two-hole am- 
plitude for holes of opposite spins in the degenerate (x, y) 
orbitals. The only contributions are the second-order di- 
agrams of Figure 1. 



the cluster results from the other. With this preparation, 
we shall be ready to present in Sect.V the general canon- 
ical transformation leading to the effective Hamiltonian 
for pairs po| . This is free from the limitations of pertur- 
bation theory and differs from the results of Sect. IV in 
an important way, although it is seen to reduce to them 
in the appropriate limit. It will turn out that the effec- 
tive interactions are additive, and the theory is able to 
accommodate the effects of phonons, off-site interactions 
interplanar coupling |21| and whatever physics one 
wants to include. The final result is an integral equation 
for pairs in the full plane, which is suitable for a numeri- 
cal study. Sect. VI is devoted to the method we used to 
solve the integral equation in finite supercells and than 
taking the asymptotic limit. Two kinds of bound states 
of different symmetries result, and Sect. VII presents 
the dependence of A on the filling and other parameters. 
Finally, Sect. VIII summarizes our conclusions. 



II. BAND STRUCTURE 



b) 



FIG. 1. The second order diagrams for the Two-Hole Ampli- 
tude. 



The system makes virtual transitions to 4-body (3 
holes and 1 electron) states. We demonstrated that this 
produced an effective interaction, which pushes down the 
singlet and up the triplet by |A( 2 )|. In this way, A can 
be redefined without any reference to the ground state 
of clusters with a different number of holes, and we are 
free from the objections based on the Jahn- Teller dis- 
tortion of odd- TV clusters. Good agreement between the 
second-order calculation and the numerical exact diag- 
onalisation results supported the interpretation. Thus, 
the cluster calculations |l6|-|lg|] showed that W=0 pairs 
are the "bare" quasiparticles that, when "dressed" ', be- 
come a bound state. That approach is inherently limited 
by the small size of solvable clusters, but allows a very 
explicit display of paired hole properties, that even show 
superconducting flux-quantization [l5|Jl9[| . Our second 
task here is to extend the theory to the full plane, fol- 
lowing the hint provided by the diagrams. In Sect. II we 
review the one-body properties of the model, in order to 
estabilish some notation and several useful formulas. In 
Sect. Ill we show how W=0 pairs in the full plane can be 
obtained for any filling. The effective interaction between 
holes is calculated by second-order perturbation theory 
in Sect. IV, with the aim of showing the relation of the 
present formalism to Cooper theory from one side and to 



The eigenstates of 

iT Q V (l/) (fc,r) = e (v) (k)ip (v) (fc,r) 



(6) 



where v is a band index, are written according to Bloch's 
theorem 



^ M 0,r) = e lkr (t> (,j) (fc,r), 



(7) 



where <fi is periodic. If r = 0, r = a and r = b are the 
positions of Cu and of the two Oxygens, respectively, the 
non-bonding band is characterized by <p^ nb ^> (k, 0) = 0; 
one then finds 



^ (nb) (M) 



cos 



N c (cos 2 + cos 2 



(8) 



^ (k,b) 



cos 
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jNc (cos 2 [% 


i] + 


COS 2 


2 


) 



(9) 



where Nc is the number of cells in the crystal and d the 
lattice parameter. For the other bands, one obtains the 
eigenvalue equation 



(±) (Jfe) - ^ 



e n ±r 



(10) 



(+ for the antibonding and - for the bonding band), with 

cos- — ros- — (11) 



r= W£ 2 + 16t 2 



and the Cu amplitudes 
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(±) 



(M) 



r(±) 



N C (2 £ (±) - e p ) ' 
while the O amplitudes are obtained from 

(e P -e M (*)) H (M) 

H (M) = o, 



(12) 



-2t cos 



kyd 



(13) 



and 



(e P -e W (fc)) 0<">(M) 
2icos f^l <f>^ (fc,0) = 0. 



(14) 



Setting x = k x d, y — k y d, we get the contours of con- 
stant energy e in the Brillouin Zone (BZ): 



cos (x) + cos (y) = q(e) . 



2. 



(15) 



(16) 



In the bonding band, q = 2 at the bottom and q = —2 
at the top; at q = the contour is a square with vertices 
at 2* (±1,0) and 2*(0,±1), where the Fermi velocity 
vanishes. Let D be the domain defined by e < Ep in the 
BZ. 

Below half filling , 



/•- o(0) r-a{x) 

dxdy = dx dy, (17) 

L) Ja{0) Ja(x) 



where a(x) = &rccos(q(Ep) — cos(a^)); at half filling 



7T pTZ— \x\ 

dxdy — I dx dy, (18) 

— 7T J— TT + \x\ 



D 



and above half filling, putting b(x) — &rccos(q(Ep) 
cos(x)) , we find 



where 



and 



dxdy = I\ + h + h , 



MO) ,7T 

h = dx dy, 

J-b(0) J-tt 



/■- 6(0) ra(x) 

h = dx dy, 

J — 7r J —a(x) 



h = / rfa; / rfy, 

/6(0) J-a(x) 



(19) 



(20) 



(21) 



(22) 



The number of holes per spin per unit cell below Ep 
in the bonding band above half filling is 



N(E f 



2vr 



d 2 k6{E F - e{k)) (23) 



differentiating, one finds the density of states 
, . e„ — 2e dx 



2ttH 2 



3 [« +1 l yd - (g- cos (x)f 



= . (24) 



where q = q(s) . Taking care of the singularity of the 
integrand at the lower limit, this expression is convenient 
for the numerical evaluation, and shows a logarithmic 
singularity at half filling (Van Hove singularity). 



III. W=0 PAIRS 

We change representation with 

cUr i ) = J2i,^ > (k,r i ), (25) 

where k runs over the BZ and v over bands. Thus, 

W = J2 4c uVl ,+ C k2,U2,+c{ 3>U3i _C kirl/ii - 



*U^~^ {kxMMM) , (26) 



Using a shorthand notation, where the band indices are 
understood, 



U (m,n,p, q) = ^2' 



— i(m-\-n—p—q)ri 



x (f>* (to, i) (f>* (n, i) (/> (p, i) (j> (q, i) 



(27) 



is the lattice Fourier transform of a periodic function; the 
momenta are limited to the BZ, so, introducing reciprocal 
lattice vectors G we may write 

U (to, n, p, q) — Nc S (to + n — p — q — G) 



G 



xU c (m,n,p,q;G), (28) 



where U involves a sum over a single cell: 

Cell 



U c (m,n,p,q;G) = Up 



-iGr t 



x 4> (to, i) 4> (n, i) 4> (p, i) <j) (q, i) , 



(29) 



where stars have been dropped since (f/s are real. The 
umklapp processes involving the basis vectors of the re- 
ciprocal lattice produce a - sign on one of the Oxygens. 
The matrix elements ( p8|) are proportional to iV^ 1 , like 
any bona fide two-body interaction. 

Still omitting the band indices, we shall mean 
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d[k] 



-k. 



J J 



\vac > 



(30) 



to be a two-hole determinant al state derived from the 
k eigenfunctions. Since k,r) — 0(fc,r), which is 
required by time-reversal symmetry, the combination 
d[k] + d[—k] is singlet and d[k] — d[— k] is triplet. 

The point symmetry Group of the Cu-0 plane is Ci v , 
and its character Table is shown in Table I. 

We introduce the determinants Rd[k] — d[Rk] = 
d[kn], R £ Ci v , and the projected states 



Triplets of the other Irreps are obtained from deter- 
minants of the form ||fci,fc 2 ||, but they vanish when 
fc 2 = — k\. If O is identihed with W, the determinantal 
matrix element does not depend on the sign of the com- 
ponents of k' , because the <fi ' s depend on k only through a 
cosine. The only thing that matters is that C4, C^,a[ and 
<j' 2 exchange k' x and k' y , while the other operations leave 
them in place. Consequently, We — Wc 2 = W G:c = W a 



and Wc 4 = W c 



W„i , and 



(31) 



W( 1 A 2 ) = W( 1 B 2 ) =0. 



(38) 



where x {R) is the character of the operation R in 
the Irreducible Representation (Irrep) 77. In the non- 
degenerate Irreps, the operations that produce oppo- 
site kfi have the same character, and the correspond- 
ing projections lead to singlets. Let Ri,i = 1, ..8 denote 
the operations of C^ v and fc, k' any two points in the 
BZ. Consider any two-body operator O, which is sym- 
metric (R,l<DRi = O), and the matrix with elements 
O i;j =< d[k]\RlORj\d[k'] >, where k and k 1 may be 
taken to be in the same or in different bands. This ma- 
trix can be diagonalized by Group Theory. Indeed, for 
each Irrep 77, consider the normalized vector with com- 
ponents ~jgX (Ri) j this is an eigenvector of the op- 
?t 



These are W=0 pairs, like those studied previously [L5|, 
but with an important change, since in the cluster calcu- 
lations the symmetries of W=0 pairs were found |l9| to 
be 1 i?2 and 1 A\. The reason for this change is a twofold 
size effect. On one hand, 1 A\ pairs have the W=0 prop- 
erty only in the small clusters, having the topology of a 
cross, and belonging to the S4 Group, but do not gen- 
eralize as such to the full plane, when the symmetry is 
lowered to C^ v ; on the other hand, the small clusters 
admit no solutions of 1 A2 symmetry at all. 

One necessary condition for pairing in clusters is that 
the least bound holes form such a pair, and this dictates 
conditions on the occupation number. In the full plane, 
however, W=0 pairs exist at the Fermi level for any fill- 
ing. The explicit W=0 singlet pair states are: 



erator matrix i^OR , as we can check by noting that $ » , = X» { [fc( _ }] (fc } ^ (fe } 

hE i , i X^(Ri)Ri 1 OR j X W (R j ) = OP* (a) 5 (a, P); ' 



hence, R\ORj is diagonal on the basis of the symmetry 
projected states. The square of the projection operator is 
P 2 (a) = V8-P (a) . Now, taking the matrix element be- 
tween the k and k' determinants, we get the eigenvalues 
in terms of determinantal matrix elements: 



V2 l 

- cos [k Ci (n - r 2 )} (k Ci , ri) <p (k C4 ,r 2 ) 
~ cos [k a (ri - r 2 )] <i> (k a , n) (f> (k a , r 2 ) 
+ cos[fc (T /(ri -r 2 )]</)(fc (T /,ri)(/)(fc (J /,r2)} (39) 



and 



0( v ,k,k') = J2x (n) (R) o R {k.k>) 



(32) $ lj42 [k, n ,r 2 ] = ^{cos [k(n - r 2 )] cf> (k, n) (fc, r 2 ) 



where 



R (k,k') = (d[k]\6\Rd{k'\). 
Thus, omitting the fc, fc' arguments, we get 
O ( 1 A 1 ) = E + C2 + Ci + O c s 



-O am + CL + 



O ( 1 A 2 ) = E + C2 + Ci + O c} 

-o„ m - o av - o a , - o a . 



O ( l B 2 ) =O e + Oc 2 

-Oa. - O a .. 



o Ci - o c} 

- Oai + 0„> 



(33) 



(34) 



(35) 



(36) 



and for the x component of the degenerate triplet 

O ( 3 E X ) =0 E - Oc 2 + 0^ - O ay . (37) 



V2 

+ cos [kciin - r 2 )}(f>(k Ci ,r 1 )(l)(kc 4 ,r2) 
-cos [k a (ri - r 2 )]0(fc o -,ri)<^(fc o .,r 2 ) 
-cos[fc (T /(ri - r 2 )](j){k <jl ,ri)4>(k a ,,r2)}, (40) 

where xo is a singlet spin function. Using the above ex- 
plicit Bloch states we can verify that both vanish for 
n = T2- This is another way to see that they have 
the W=0 property. All the distinct W=0 pairs of each 
symmetry may be labeled from the fc points in 1/8 
of the empty part of the BZ: for instance, those with 
k x > k y > 0. We shall denote such a set by e/8. 



IV. CANONICAL TRANSFORMATION-SECOND 
ORDER 

In this Section, we shall parallel as closely as possi- 
ble the Cooper theory [g2| where an effective interac- 
tion involving phonons is introduced via an approximate 
canonical transformation. The treatment we give here 
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is based on second-order perturbation theory and is not 
yet adequate for the present purposes, but is useful as a 
prelude to the more complete theory given in the next 
Section. We wish to stress the formal analogy with the 
BCS theory and make contact with the A expression de- 
rived from cluster studies. Suppose we add two holes to 
a background Fermi sphere. We use Roman indices for 
the zeroth-order pair states, which satisfy 



H \m) = E v 



(41) 



with E m — 2e (m) . The basic idea here is that of obtain- 
ing an effective interaction for two holes using the process 
of the diagrams in Fig.l. Accordingly, the pair states are 
coupled to the set of 3 hole-1 electron intermediate states, 
which satisfies 

H \a) = E a \a) . (42) 
Thus, the effect of the perturbation is of the form 

W\m) = Y,\m')W m ,, m + Y.\ a ) W °<<™ ( 43 ) 

m' a 

where W m ' tm stands for the interaction between pairs, 
while a is the set of 3 hole-1 electron intermediate states 
which are coupled to the pairs by the diagram; 



W\a) 



E 



m) W„ 



(44) 



this is accomplished if 

W a . m + (E a - E m ) A 







(48) 



with all others A-matrix elements equal to zero. 

Assuming that the denominators do not vanish (more 
about that later) we obtain 



(a\ A\m) = 



(a\ W \m) 

E m — E n 



(49) 



In H, the holes interact through the effective vertex of 
Figure 2. 



a) 



The off-diagonal W m ^ m elements vanish for W=0 pairs; 
on the other hand, in the many-body problem, the di- 
agonal elements W m ^ m do not vanish, because of the 
tadpole-diagram contributions to the self-energy of the 
holes. However, the diagonal terms do not contribute 
to the effective interaction, but simply renormalize the 
e parameters in H - We had already met these tadpole 
diagrams in the calculation of A^ 2 ) for clusters where 
they were seen to cancel in the calculation of (H). Any- 
way, we keep W m i >m in (^3|) which allows to introduce 
the effects of other interactions which are not included in 
the Hubbard model. In the present Section we shall con- 
sider IF as a small operator and look for an approximate 
canonical transformation such that the new Hamiltonian 
H decouples the a states to first order; then, H operates 
on the space of pairs. In other terms, we seek a first-order 
anti-Hermitean operator A such that 



H = e- A He A 



(45) 



has no linear term connecting the |a) states with the \m) 
states. Since 

H = H+[H,A] + ^ [[H, A] , A] + (A 3 ) (46) 
expands to read 



H = H+ [H , A] + [IF, A] + - [[H , A] , A] + ... (47) 



FIG. 2. The effective hole-hole vertex (blob) defined by the 
second-oder H. The dashed lines represent W interactions. Linked 
diagrams belong to W e ff, and unlinked ones to F. 

The interaction lines in the right hand side represent 
IF; while a) and its specular counterpart b) correspond 
to the effective interaction diagrams of Figure 1, c) and d) 
are unlinked and represent self-energy processes of each 
hole. The same situation occurs in Cooper theory, where 
the self-energy processes involve the emission and absorp- 
tion of a phonon by the same electron. 

We may write 



H = H a + F + W eff 



(50) 



where F is diagonal in the pair space, like H , and cor- 
responds to the unlinked self-energy diagrams, while the 
effective interaction operator is W e tf. Like in Cooper 
theory, F will be dropped. Formally, we write 



F + W eff 



and 



2(m\F + W eff | to' 
1 1 



(51) 



(52) 



G 



The a states are 3 hole-1 electron determinants which 
carry no quasi-momentum. We write 

\a >= | \\(k' + q + k 2 ) + ,k 2 -,-q-,-k'_\\ > (53) 

where k 2+ is the electron state and pedices refer to the 
spin direction; a states with opposite spin indices arise 
from the specular diagram, contribute similarly and yield 
a factor of 2 at the end. The unperturbed eigenvalues are: 



E n 



s(k' + q + k 2 )-e (fc 2 ) + S (q) + E (k') . (54) 



Let us first do the calculation on a determinantal basis. 
We work out the interaction matrix element: 



< \\{k' + q + k 2 ) + ,k 2 -,-q-,-k'_\\ \W\d[s] >= 
S(q- s) U{q + k' + k 2 , -k', s, k 2 ) 
-5 (k' -s)U(q + k' + k 2 , -q, s, k 2 ) 



(55) 



two of the 4 terms in the product in (|52|) are propor- 
tional to S (p — s) and belong to F, while the cross terms 
yield identical contributions to the effective interaction. 
Symbolically, we may write 



(ml W, 



■ J f l m '> = 2 H [ W m,aW a , m '} Linked 



E,, 



E„ 



(56) 



Including the specular diagram by a factor of 2, we write 
the interaction between determinantal states (with s and 
p empty). 

occ 

(d[p}\W eff \d[ S })^2j2e(e{s+p + k)-E F ) 



xU(s 



p + k 
{ 



-p, s, k)U {p, fc, s + p + fc, — s) 
1 



e (s + p + k) — e (k) — e (s) + e {p) 

I }. 

e (s + p + k) — e (k) — e (p) + e (s) 



(57) 



The interaction between symmetry projected states is ob- 
tained using (|3^j33]) and is given by 



(® a \p]\W eff \$ a [s]} 



J2x {a) (R)(d[ P ] 



w, 



eff 



Rd\ 



(58) 



The above theory applies not only to the plane but also 
to clusters, provided that they are fully symmetric and 
allow W=0 pairs. The natural basis set is provided by 
the orbitals, rather than the Bloch states; in the small 
clusters with open boundary conditions that we studied 
previously, pairing occurs when a degenerate state of the 
Irrep e(x,y) is occupied by two holes. Then, in terms of 
orbitals, the only 1 B 2 W—0 pair available for pairing is 



$1R, = 



\\x+y-\\ + \\y+x-\\ 

V2 



(59) 



The diagonal element of W e ff on this state gives the 
second-order energy shift of the pair; one calculates the 
W matrix on the orbital basis and after some algebra one 
finds the above A^ 2 ) expression (^). In the small clusters, 
the diagonal element provides the binding energy, while 
the localization is granted by the boundary conditions. In 
the full plane, there is a continuous infinity of W=0 pairs 
for each of the allowed symmetries, and the off-diagonal 
elements are essential to form a localized wave packet. 
Therefore the localized state is not known a priori, and 
A cannot be obtained as an expectation value of W e f / , 
but must be found by solving a two-body problem with 
interaction W e ff (See Sect. VI). 

In the continuum limit, the integral in ( p38| ) must be un- 
derstood as a principal part, because the denominators 
vanish in a set of vanishing measure within the integra- 
tion domain. While in principle this is no trouble, the 
numerical evaluation is difficult. However, this compu- 
tational difficulty is an artifact of the perturbation ap- 
proach, as will become clear in the next Section. 

Above all, the approximate canonical transformation is 
open to question because it assumes that W be weak. In 
small clusters, it works remarkably well, even for at large 
couplings, but one could argue that the reason might be 
that A is still small compared to the average distance 
between the discrete levels. When the levels form a con- 
tinuum, the applicability of perturbation theory is not 
obvious. Therefore, it is important to put the theory on 
a clearer and firmer basis, as we do in the next Section. 
A posteriori, it turns out that high order diagrams mer- 
erly renormalize some quantities without destroying the 
structure of the theory, because it is a property of W=0 
pairs that the interaction is dynamically small. 



V. CANONICAL TRANSFORMATION 
ORDERS 



ALL 



Suppose the Cu-0 plane is in its ground state with 
Fermi energy Ep and a couple of extra holes are added. 
We wish to show that by a canonical transformation |2(J 
one obtains an effective Hamiltonian which describes the 
propagation of a pair of dressed holes, and includes all 
many-body effects. 

The exact many-body ground state with two added 
holes may be expanded in terms of excitations over the 
vacuum (the non-interacting Fermi sphere) by a configu- 
ration interaction: 



(60) 



here m runs over pair states, a over 4-body states (2 holes 
and 1 e-h pair), (3 over 6-body ones (2 holes and 2 e-h 
pairs), and so on. To set up the Schrodinger equation, 
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we consider the effects of the operators on the terms of 
|*o >• We write: 

H \m >= E m \m >, H \a >= E a \a >, ... (61) 

and since W can create or destroy up to 2 e-h pairs, 

W\m >= Y^ w m', m \m' > +El a > W a , m 

m' a. 

+ ]T|/3 > W 0,m- (62) 



As explained in the previous Section, W rn ^ m does not 
contribute to the effective interaction for W=0 pairs in 
our model; however we keep it for generality, since it al- 
lows to introduce the effect of the exchange of phonons 
and other quasiparticles that we are not considering. For 
clarity let us first write the equations that include explic- 
itly up to 6-body states; then we have 

W\a >= > W m , a + J2\ a ' > W a ', a 

m a' 

+£2\p > Wp, a (63) 


where scattering between 4-body states is allowed by the 
second term, and 

w\p >= W) w m >, fi + e i«> Wc^p 



J2 W) 

0' 



(64) 



In principle, the Wp> p term can be eliminated by taking 
linear combinations of the complete set of (3 states: when 
this is done, we get a self-energy correction to Ep and a 
mixing of the vertices, without altering the structure of 
the equations. The Schrodinger equation yields equations 
for the coefficients a, b and c 

(E m — Eq) 

+Y am ' Wm - m ' + Y haWm > a + Y c Wm >0 = ( 65 ) 



(E m — Eq) a m + y^a T , 

m ' 

+E 5 « 



w„ 



E 



Eq — Ep 



W„ 



E 



W m ,pWp, a 
Eq — Ep 



= (68) 



(E a - Eq) b a + Y a ™' Wa 



-5> 



w ataf + J2 







Wp,g>W a ,p 

Eq — Ep 



= 



(69) 



Thus we see that the role of 6-body states is just to renor- 
malize the interaction between 2-body and 4-body ones, 
and for the rest they may be forgotten about. If Eq is out- 
side the continuum of excitations, as we shall show below, 
the corrections are finite, and experience with clusters 
suggests that they are small. Had we included 8-body 
excitations, we could have eliminated them by solving 
the system for their coefficients and substituting, thus 
reducing to the above problem with further renormaliza- 
tions. In principle, the method applies to all the higher 
order interactions, and we can recast our problem as if 
only 2 and 4-body states existed. Again, the W a ', a term 
can be eliminated by taking linear combinations of the 
a states: when this is done, we get a self-energy cor- 
rection to E a and a mixing of the W m , a vertices. The 
equations become 

(E m - E a ) a m + Y J a m'W m , m ' + Y haWm > a = ( 70 ) 



(E a - E ) b a + Y a m'W atm ' = 



(71) 



Solving for b a and substituting in the first equation we 
exactly decouple the 4-body states as well. The eigen- 
value problem is now 

(E - Em) a m = E a ™' { w m,m' + (m\S[E ]\m)} , 



(72) 



(E a — E )b a 

+Y, a m>W a , m i + J^fta'Wa.o' + ^CpW a ,fi = (66) 
m' a' (3 

{Ep - Eq) Cp + J2 a m'Wp, m ' + J2 b <*' W 0,<*' = ( 67 ) 



where Eq is the ground state energy. Then, we exactly 
decouple the 6-body states by solving the equation for cp 
and substituting into (B3,pq), getting: 



where 



i5[^]K)=E <m|wl rv wK> - <™ 



Equation (|72j) is of the form of a Schrodinger equation 
with eigenvalue Eq for pairs with an effective interac- 
tion W + S. Then we interpret a m as the wave function 
of the dressed pair, which is acted upon by an effective 
Hamiltonian H. The change from the full many-body 
H to H is the canonical transformation which general- 
izes the one of the previous Section to all orders. In the 
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new picture, the holes interact through an effective ver- 
tex with infinitely many contributions, some of which are 
shown in Figure 3. The linked contributions represent 
repeated exchange of electron-hole pairs, and may con- 
tain self-energy insertions; all these contributions make 
up the effective interaction. The unlinked diagrams are 
pure self-energy. Thus, the scattering operator S is of the 
form S = W e ff + F, where W e f / is the effective inter- 
action between dressed holes, while F is a forward scat- 
tering operator, diagonal in the pair indices m ,m' which 
accounts for the self-energy corrections of the one-body 
propagators: it is evident from ( |72| ) that it just redefines 
the dispersion law E m , and, essentially, renormalizes the 
chemical potential. Therefore F must be dropped, as in 
Cooper theory and in the last Section. 



:d 



+ q: 



FIG. 3. Some terms of the expansion of the effective hole-hole 
vertex (blob) defined by H. The dashed lines represent W inter- 
actions. Linked diagrams belong to W e ff, and unlinked ones to 
F. 



So, the effective Schrodinger equation for the pair reads 



(H Q + W + W eff )\a>= E Q \a> 



(74) 



and we are interested in the possibility that E = 2Ep — 
|A|, with a positive binding energy |A| of the pair. Any 
other pairing mechanism not considered here, like off-site 
interactions, inter-planar coupling and phonons, can be 
included as an extra contribution to W m > >m which just 
adds to W e ff. 

We emphasized the fact that the canonical transforma- 
tion is exact because in this way our argument does not 
require U/t to be small. In the numerical calculations 
below, we shall drop W and approximate W e ff neglect- 
ing the 6-body and higher excitations; in other terms, 
we shall compute the bare quantities. At least the struc- 
ture of the solution is exact when expressed in terms 
of renormalized matrix elements, and the process can be 



systematically improved. The a states are those of Equa- 
tion ([53J) and the interaction matrix element is given in 
Equation (p5|). Working out (|7^) we find that the prod- 
uct in the numerator yields 4 terms; two are proportional 
to S(p — s) and belong to F, while the cross terms yield 
identical contributions to W e f / . Using ( ^2|3^ ) , we obtain 
the effective interaction between W=0 pairs: 

(*„ Ml w e// 1$„[«]> = 

occ 

4 (r) ( e ( Rs +p+ k )- E F) 

U(Rs+p + k. -p, Rs, k) U (p, k.Rs + p + k, -Rs) 



e(Rs+p + k)-e(k)+e (s) + e (p) - E 



(75) 



The sum is over occupied k with empty Rs + p + k. 
Note that W e // does not depend on the sign of U. We 
may see that the perturbation theory of Section V yields 
the arithmetical mean of the two unperturbed limits 
Eq — > 2e (p) and Eq — > 2e (s) . The diagonal elements 
of (|75| ) are also clearly related to the A^ 2 -* expression de- 
rived from perturbation theory for the fully symmetric 
clusters which are special cases of the present the- 
ory. The expressions (|o]) , (|58|) and ( |75|) are characterized 
by the symmetry-induced quantum mechanical interfer- 
ence of several terms. The sum can be positive or nega- 
tive depending on the Irrep. This interference produces a 
partial cancellation, and the absolute value of the result 
is typically much smaller than individual contributions. 
This explains why the interaction is dynamically small 
for W=0 pairs: they have no direct interaction, and be- 
cause of the interference the effective interaction is re- 
duced compared to what one could expect by a rough 
order-of-magnitude estimate. However, the presence of 
the 9 functions and the anisotropy of the integrands pre- 
vent a total cancellation. 

In (|^) , the vanishing of the denominators is prevented 
by the discrete spectrum of the cluster; however, when 
applying ( |58|) to the plane the integrand is singular. This 
complication disappears in the full expression ( f75| ) since 
we are interested in bound pairs; then, Eq is below the 
continuum and the denominators never vanish. 

The s and p indices run over 1/8 of the BZ. We denote 
such a set of empty states e/8, and cast the result in the 
form of a (Cooper-like) Schrodinger equation 

e/8 

2e (fc) a (k) + J2 w eff {k, k') a(k')=E a (k) (76) 
ft' 

for a self-consistent calculation of Eq (since W e /f de- 
pends on the solution). Let Nc be the number of cells 
in the crystal. The U matrix elements (f28|,[29|) scale as 



Nq and therefore W e ff scales in the same way. For an 
infinite system, Nc — > 00 , this is a well defined integral 
equation. The existence of a clear-cut continuum limit 
means that the problem is well posed, but this equation 
is very difficult to treat analytically. The interaction ma- 
trix is complicated by the Umklapp discontinuities, other 







discontinuities come from the limitations to occupied or 
empty states, and there are several independent variables 
since everything is anisotropic. Therefore, we must re- 
sort to a numerical treatment. For the sake of simplicity, 
we shall neglect the minor contributions from the higher 
bands and consider the dominant intra-band processes, 
in which empty states belong to the bonding band. 



VI. SUPERCELL CALCULATIONS 

Several good supercell calculations devoted to the 
problem of pairing have been reported to date |23f| , but no 
conclusive evidence was reached, because of the difficulty 
of dealing with size effects. Any brute force approach will 
miss the result by far, because size effects are quite im- 
portant, or if you wish, the minimum size useful for simu- 
lating the plane is enormous. Now we have the advantage 
of starting from an analytic result (Equation (|75|)) which 
holds for the full plane. However the analytic solution 
is difficult, so we must discretise the problem. The best 
way is to work with a supercell of Nsc x ^Ysc = Nc 
cells, with periodic boundary conditions. The advantage 
remains, however, because Equation (|75|), however hard, 
is easy compared to direct diagonalisation. We are able 
to solve the problem in a virtually exact way for Nc = 30 
or 40, while the calculation becomes somewhat hard with 
larger supercells. 

We loooked for triplet bound states, but we did not find 
them, because, as in the clusters, W e ff is repulsive for 
triplets. Therefore the rest of the analysis is concerned 
with the W=0 singlets, which do show pairing, as we 
demonstrate below. 

The A < behaviour that we obtain in supercell cal- 
culations is_in line with our previous findings in small 
©, @, @, " 



clusters ( JL5|, |6j, fl7|], (18)). The supercells have 
periodic boundary conditions instead of the open bound- 
ary conditions we were using previously and are very 
much larger systems than those that one can afford by 
exact diagonalisation. However, we are interested in the 
asymptotic behavior for Nc — > oo and A depends on U's 
and Nc in a peculiar way, within the range of attain- 
able supercell sizes. It increases with Nc for large U's 
and drops for small U's; we must understand the depen- 
dence in order to be able to make reliable extrapolations. 
Cooper p2l assumed a constant negative interaction — Vc 
within a distance ljd from Ep, where ljd is the Debye 
frequency. His result is well known: 



|A r 



2lo d 



e v cPF — 1 



(77) 



where pp is the density of states at the Fermi level. 
The orders of magnitude are u>d ~ 25 meV,pF « 
0.3eV _1 /cell and Vc ~ leV/cell, and |Ac| is in the meV 
range. Here, no LUp> is involved, and all the empty states 
contribute; however, we recall that all the distinct W = 
pairs may be labeled from the k points in the e/8 set. We 



define the Uniform Interaction Model (UIM) in which a 
constant negative interaction — V prevails for k and k' in 
e/8. Both the Cooper and the UIM models can be solved 
in supercell calculations, and show the same qualitative 
behavior with increasing Nc- For small supercells, the 
binding energy is of the same order as V, while for large 



cells it tends to the asymptotic value — A 



sympt 



which 



may be very different. This is similar to the behavior we 
observe in the Hubbard model. Moreover, we can com- 
pute A for a given filling and Nsc according to (^6|) , and 
determine the value of V which gives the same A in the 
UIM. This will be V e ff, that is, the effective V of our 
theory, characterizing the strength of the attraction by a 
single quantity. Fortunately, V e ff does not depend too 
much on the supercell size, and is seen to converge with 
increasing Nsc- 

We use as input data the current estimates (in eV) 
t=1.3, e p =3.5, £<j=0, U p = 6s, Ud — 5.3s, where s is 
a scale factor induced by renormalization. At s=2.121, 
with Ep=-1.3 eV, we get for 1 B 2 pairs the results shown 
in Tabic II. Here, n to t is the filling and V e ff is derived by 
comparison with UIM calculations. We see that although 
A | decreases monotonically with increasing supercell size 
up to Nsc = 40, V e ff is fairly stable and is clearly not 
dropping to zero; it corresponds to |A asympt | values of 
about 20 meV. 

In addition, the UIM can be solved in the thermo- 
dynamic limit; writing Eq = 2Ep — |A| we can esti- 
mate A. asymp t for Nc — * oo . Indeed, the Cooper-like 
Schrodinger equation (|7^ ) with W e ff — — jf^, V > 0, 
leads to 



V 



dep (e) 



2(e-E F ) + \A a 



sympt | 



(78) 



where p is the density of states. This is readily solved 
numerically. The dependence of A on V is qualitatively 
similar to the dependence of Ac on Vc, and good fits 
can be obtained using the Cooper formula with effective 
Lop) and pp parameters, like those shown in Table III. 



VII. DOPING DEPENDENCE 

In Table IV, we consider 1 B 2 pairs at Ep = — 1.2 eV 
(with the current parameters, half filling corresponds to 
Ep = —1.384 eV, so the filling we are considering is 
n-tot ~ 1-3 at Nsc — 30). We report values of A and V e ff 
as obtained from supercell calculations at -/Vsc=12,20 
and 30 for 3 values of the scale factor s. The extrap- 
olation to the infinite supercell is obtained from the UIM 
with the Veff value obtained at Nsc = 30. We see that 
the relatively mild Nsc dependence of V e f / supports the 
use of the UIM to extrapolate the results to the thermo- 
dynamic limit, and there is a clear indication of pairing 
with sizable binding energies. The s dependence of V e tf 
is roughly linear, while A depends exponentially on s. 
Table V presents the results for 1 B 2 pairs at Ep = —1.1 
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eV ( n to t ~ 1.4 at ATgc = 30). The trend is similar, but 
V e ff is seen to increase with doping. Tables VI and VII 
show the results for the 1 A 2 pairs. These are seen to lead 
to bound states as well, with comparable A values; the 
trend with doping is opposite, however, and the gap is 
nearly closing at Ep = —l.leV. A necessary condition 
for superconducting flux quantization is that two kinds 
of pairs of similar binding energy and different symme- 
tries exist A similar conclusion was reached inde- 
pendently by other Authors juj. Moreover, evidence of 
mixed (s + id) symmetry for the pairing state has been 
amply reported in angle-resolved photoemission studies 
p4j . This remark leads to the prediction that in this 
model superconductive pairing disappears with increas- 
ing ntot, while a different sort of pairing prevails; in real- 
ity the Cu-0 plane prefers to distort at excessive doping, 
and in a distorted plane the present mechanism, based 
on symmetry, could be destroyed. 

VIII. CONCLUSIONS 

We propose a general framework for the effective inter- 
action between two holes, based on the three-band Hub- 
bard model but ready to include extra interactions as 
well. An effective Hamiltonian can in principle be ob- 
tained by a systematic canonical transformation includ- 
ing any kind of virtual intermediate states. We obtained 
the closed-form analytic expression of the effective inter- 
action including 4-body virtual states. This describes re- 
peated exchange of an electron-hole pair. The argument 
does not depend on perturbation theory, and the equa- 
tions retain their form, with renormalized parameters, at 
all orders. The previous exact-diagonalisation results of 
cluster calculations are special cases. The resulting in- 
tegral Equation ([76|), with the effective interaction(f75|), 
is valid for the full plane. Since an analytic treatment is 
prohibitive, we resort to a numerical treatment by super- 
cell calculations. Even so, the solution is hard. Depend- 
ing on the parameters, extremely large supercells may be 
needed to obtain the convergence of the pair energy A to 
the bulk limit; however, we find that precisely the same 
effect occurs in the Uniform Interaction Model where a 
constant effective interaction V is assumed. We define 
V e ff as the value of V that inserted in the UIM yields 
the same A in supercell calculations as our integral equa- 
tion. Since V e f / converges to the bulk value much more 
readily, we are able to go to the asymptotic limit, and to 
show the instability of the Fermi liquid in the model at 
hand. A asymp t values in the range of several tens to a 
few hundreds of meV are obtained if we multiply the U 
parameters by a scale factor s which is somewhat larger 
than 1. The values Ud — 5.3 eV, U p = 6. eV differ ap- 
preciably from other literature estimates [p5[ , and must 
depend on the compound and doping. For example, in 
La2CuC>4, Up=4 eV and £7^=10.5 eV have been recom- 
mended p6j. However, since the screening excitations are 



explicitly accounted for in the Hamiltonian, it is reason- 
able that the input U's must be somewhat larger than 
the fully screened interaction. Moreover, contributions 
from phonons and other mechanisms can be included as 
additive terms of W, and must be relevant for a compar- 
ison with experiment. We find that 1 A-i pairs are more 
tightly bound close to half filling, but 1 £?2 pairs are fa- 
vored when the filling increases. We remind the reader 
here that these symmetry labels are not absolute, but 
depend on the choice of a gauge convention |^7j . We get 
attraction and pairing at all fillings, but the binding en- 
ergy of the l A2 pairs drops by orders of magnitude as 
the filling increases; thus, there is no chance of supercon- 
ducting flux quantization too far from half filling. So, 
we do not predict that superconductivity occurs outside 
some range of hole concentration. However, pairing is 
still there, even for large doping: since the present mech- 
anism is driven by symmetry it works unless the system 
distorts. In fact, at excessive dopings the real supercon- 
ductors develop stripes, and become normal metals. The 
three-band Hubbard model might be too idealized to al- 
low a detailed comparison with experiments; however we 
stress that the approach presented is far more general 
than the model we are using, and can be applied to more 
realistic Hamiltonians. This is the main result of the 
present paper. 
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TABLE I. The Character Table of the C iv Group 



Nsc 


ntot 


—A (meV) 


Veff (eV) 


Aasympt (meV) 


18 


1.13 


121.9 


7.8 


41.6 


20 


1.16 


42.2 


5. 


9.0 


24 


1.14 


59.7 


7. 


28.9 


30 


1.14 


56. 


5.7 


13.2 


40 


1.16 


30.5 


6.6 


23.4 



TABLE II. Binding Energy of 1 B 2 Pairs in supercells. 



E F {eV) 

-1.35 

-1.3 

-1.2 

-1.1 



u D (eV) 
0.4545 
0.492 
0.520 
0.5237 



PF(eV~ 1 /Cell) 

0.0486 

0.0415 

0.0338 

0.02905 



TABLE III. Effective Cooper parameters for fitting the V 
dependence of A at several Ep values. 





s — 


V2 


s — 


3 

V2 


s — 


2\/2 


Nsc 


-A(meV) 


Veff 


-A(meV) 


Veff 


-A (meV) 


Veff 


12 


126 


6.75 


251. 


10. 


469. 


14.23 


20 


63. 


6.9 


135. 


10.6 


320.5 


15.9 


30 


27. 


6.3 


198. 


13.8 


460. 


20.8 


oo 


8.8 


6.3 


131. 


13.8 


333. 


20.8 



TABLE IV. Data for 1 B 2 pairs at E F = -1.2 eV; Veff is 
in eV. 





s — 


V2 


s — 


3 

vl 


s = 


2^2 


Nsc 


-A (meV) 


Veff 


-A(meV) 


Veff 


-A (meV) 


Veff 


12 


163. 


9. 


409. 


15.7 


688. 


22.1 


20 


151. 


10.8 


371. 


17.6 


635.5 


24.3 


30 


80. 


10.3 


270. 


17.66 


520. 


24.9 


oo 


35.6 


10.3 


170.8 


17.66 


362. 


24.9 



TABLE V. Data for 1 B 2 pairs at E F = -1.1 eV; V eff is in 
eV. 
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s — 


V2 


s — 


3 

V2 


s 


2^2 


Nsc 


A(meV) 


Veff 


A (meV) 


Veff 


\ (meV) 


Veff 


12 


141. 


7.2 


274. 


10.43 


423. 


13.4 


20 


46.2 


6.5 


96. 


9.1 


158.2 


11.4 


30 


25.9 


6.23 


64.8 


8.8 


131. 


11.6 


oo 


8.35 


6.23 


34. 


8.8 


81.7 


11.6 



TABLE VI. Data for 1 A 2 pairs at E F = -1.2 eV; V eff is 
in eV. 





s — 


V2 


s = 


3 

V2 


s = 


2\/2 


Nsc 


A (meV) 


Veff 


A (meV) 


Veff 


A (meV) 


K// 


12 


72. 


5.5 


144. 


8.4 


228. 


11. 


20 


10. 


2.4 


26.5 


4.38 


61. 


6.78 


30 


5.48 


2.5 


19. 


5.8 


45. 


8.16 


00 


1.5 


2.5 


2.6 


5.8 


14.4 


8.16 



TABLE VII. Data for X A 2 pairs at E F = -1.1 eV; V e// is 
in eV. 



